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Abstract 

We theoretically demonstrate the formation of multiple bi-stability regions in the temperature pattern on the interface between a 
cuprous oxide quantum well and DCM2:CA:PS organic compound. The Frenkel molecular exciton of the DCM2 is brought into 
resonance with the IS quadrupole Wannier-Mott exciton in the cuprous oxide by " solvatochromism" with CA. The resulting hybrid 
is thermalized with surrounding helium bath. This leads to strongly non-linear temperature dependence of the laser field detuning 
from the quadrupole exciton energy band which is associated with the temperature induced red shift of the Wannier exciton 
energy. Numerical up and down-scan for the detuning reveals hysteresis-like temperature distribution. The obtained multiple bi- 
stability regions are at least three orders of magnitude bigger (meV) than the experimentally observed bi-stability in bulk cuprous 
oxide (neV). The effective absorption curve exhibits highly asymmetrical behavior for the Frenkel-like (above the IS energy) and 
Wannier-like (below the IS energy) branches of the hybrid. 
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1. Introduction 

The proposed of enhancement of nonlinear optical re- 
sponse from the quantum confined Wannier-Mott excitons 
(WE) by means of resonant hybridization with organic 
Frenkel excitons (FE) was first stated in the works of Agra- 
novich et al. (1994); Ivanov and Tikhodeev (in press 2007); 
Agranovich and Bassani (2003). The dipole-dipole hybrid 
(DDH) is an appropriate coherent linear combination of 
large radius as 2 dipole allowed WE and small radius but 
big oscillator strength f F dipole allowed FE. The new hy- 
brid is characterized by big oscillator strength and small 
saturation density (strong coupling regime) . The expected 
hybrid optical nonlinearities are large because the ideal 
bosonic approximation starts to break down at the density 
of the hybrid close to the saturation density. Close to the 
saturation, due to overlapping between exciton wave func- 
tions, they exhibit rather strong exchange interaction and 
space filling factor. Compared to the bare WE nonlinear 
response, the hybrid response is enhanced roughly by the 
factor of the ratio of the FE and WE oscillator strength. 



In our recent work Roslyak and Birman ( arXivxond- 
mat/0703268) we demonstrated formation of a new type 
of the hybrid. It occurs between dipole forbidden but 
quadrupole allowed IS exciton in cuprous oxide quantum 
well and FE in an organic composite called DCM2:CA:PS. 
For the details on the organic see works of Madigan and 
Bulovic (2003) and Bulovic et al. (1999). The FE formed 
on the DCM2 molecules is brought dynamically into res- 
onance with the WE in the quantum well by means of 
"solid state solvation" [Baldo et al. (2001)] mechanism. 
The quadrupole- dipole hybridization (QDH) is of the same 
order of magnitude (meV) as that of the DDH: the small 
oscillator strength of the IS quadrupole exciton [Gross 
and Karryev (1952)] is compensated by its strong spa- 
tial dispersion. This QDH exhibits the key features of big 
oscillator strength and narrow line width. 

In the present work we utilize the properties of the hybrid 
to enhance and modify nonlinear effects associated with 
the WM part of the hybrid. Due to the small radius of 
the IS* quadrupole 7 A the hybrid may achieve Bose 
condensation before it reaches the critical bulk density 3 . 
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3 we leave aside many questions about the ability for actual con- 
densation of the hybrid for our future work 
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Therefore, density dependent optical non-linearities which 
were important in the DDH case are negligibly small for the 
QDH. Hence, we focus our attention on another nonlinear 
optical phenomena generic to the cuprous oxide, namely, 
the photo-thermal bi-stability effect. 

The photo-thermal bi-stability effect [Roslyak and Bir- 
man (arXiv:cond-mat/0703650v2); Dasbach (2004)] man- 
ifest itself for the density of the pumping laser as small 
as P ex — 1 mWt. Due to 'weak' interaction with acoustic 
phonons the energy of the WE will experience a red shift 
linear with the temperature [Passler (1998); Weiser (1992); 
Varshi (1967)]. The laser heating is balanced by the sur- 
rounding helium cooling. Hence, the resulting temperature 
on the interface is described by a nonlinear equation with 
more than one solution for some values of the laser detuning 
from the resonant energy level. This results in a hysteresis- 
like pattern for the temperature and absorption. 

We have shown that the fact of hybridization with the 
organic enlarges the bi-stability region from 200 neV for 
the bulk cuprous oxide to the order of meV for the hybrid 
due due to its large oscillator strength. Another remarkable 
effect of the hybridization is that the QDH exciton induces 
multiple highly asymmetrical photo-thermal bi-stability ef- 
fects associated with two dispersion branches of the hybrid. 
The effective absorption maximum for the lower branch of 
the QDH exciton experiences a red shift of the order of 
meV and strongly depends on the pumping intensity. 

In the next section we propose a nonlinear absorption ex- 
periment and we give an appropriate mathematical model 
to observe the hybrid bi-stabilty. Our last section is devoted 
to discussion of the numerical results following from our 
model, and comparison with reported results for the bulk 
cuprous oxide [Passler (1998)]. 

2. Nonlinear absorption by the hybrid 
quadrupole-dipole exciton 

Our proposed configuration is shown in the following dia- 
gram (See Fig.l). For the nonlinear absorption experiment 
let us consider a mono-layer of cuprous oxide placed upon a 
thin film of an organic compound Wc take the layer width 
L w approximately equal to the size of a unit cell a = 4.6 A. 
The organic we propose is an 'solid state solvation' of the 
red laser dye DCM2 4 embedded in a transparent host 
of polystyrene (PS) doped with the polar small molecule 
camorphic anhydride (CA) [Bulovic et al. (1999)]. One sta- 
tionary laser (probe) is tuned around the IS" quadrupole 
transitions in the cuprous oxide hw w Eis — 2.05 eV 
(WE). 

The pumping laser pulses of intensity P ex and energy 
foju = Edcm2 allow DCM2 molecular excitations. We con- 
sider the molecular excitation as a Frenkel exciton (FE). 
Therefore it can be viewed as a 2D lattice of dipolcs [i F = 
20 D, f F = 0.45 Madigan and Bulovic (2004) placed at 



z' w L w /2. After the pumping pulse the FE experiences 
a red shift linear with the CA concentration due to non- 
resonant Forster energy transitions to the CA molecules 
(See Fig.l). 
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Fig. 1. Schematic representation of a possible experimental set-up 
to nonlinear absorption by quadrupole-dipole exciton. Here the in- 
organic Cu^O quantum well provides the 15 quadrupole WE. The 
DCM2 part of the organic "solid state solute" provides dipole al- 
lowed FE; the PS host prevents wave function overlapping between 
organic and inorganic excitons; CA under proper concentration al- 
lows tuning of the excitons into resonance. 

To achieve the red spectral shift of the FE into the res- 
onance with the quantum confined WE, one has to adjust 
the CA concentration to 22% [Roslyak and Birman ( 
arXiv:cond-mat/0703268)]. To avoid complicated problems 
of the dynamics of the hybridizaton wc assume that the FE 
and WE are in exact resonance once the DCM2 enrgy is 
in close proximity to the WE enrgy Edcm2 — Eis < IV 
Now, we introduce the quadrupole-dipole coupling param- 
eter [Roslyak and Birman ( arXiv:cond-mat/0703268)]: 
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Here ap, o-b are the Bohr radius of the FE and WE; e and 
£ are the corresponding bulk dielectric constants of the or- 
ganic and cuprous oxide. The quadrupole transition ma- 
trix element Q xz may be estimated from the corresponding 
bulk oscillator strength /[no] = 3.6 x 10~ 9 via the following 
identity [Moskalenko and Libcrman (2002)]: 
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4 [[2-mcthyl-6-2-(2,3,6,7-tetrahydro-lH, 5H - benzo[i,j] - quinolizin 
- 9 - yl) - ethenyl] - 4H - pyran - 4 - ylidenc] propane dinitritle 



Here the energy gap of cuprous oxide is denoted as E g = 
2.173 eV taken at T = 1.7 K of surrounding helium; A rep- 
resents the effect of the quantum confinement, and |ko| = 
2.62 • 10 5 cra -1 . The maximum value of the coupling pa- 
rameter is about 2 meV and for the hybridization to be 



effective it must be of the same order of magnitude as the 
dissipative width of the hybrid jh. 

We assume that the cuprous oxide has purity of 99.99% 
with the reported line- width of h-fis = 0.1 meV i.e. pico- 
second lifetime [Shen et al. (1996)] while the FE life-time 
is of the order of ns. In this case the hybrid life-time is 
dominated by its inorganic part fvy k> fvyis- 

The hybridization will strongly affect the absorption of 
the probe laser and consequently, it will also modify the 
photo-thermal bi-stability characteristics of the cuprous 
oxide. To study this effect, let us assume that the whole 
system is placed into helium bath with the temperature 
Tbath — 1-7 K. The temperature of the exciton gas T at the 
interface between the organic and inorganic is determined 
by the balance between the photo-thermal heating due to 
the exciton absorption and the cooling by the helium bath. 

We will assume that the exciton gas is in equilibrium 
with the helium bath (thermalization) 5 . For thermalizcd 
excitons one can treat the temperature T as just another 
parameter of the system. 

The rate of helium cooling is given by the following linear 
term: —C(T — Tb a th)', the helium thcrmo-conductivity is 
denoted as: C = 2.8 x 10 3 s" 1 [Dasbach (2004)]. To find the 
heating term we derive below a semi-classical expression 
for the photo-thermal absorption of the hybrid. 

The Hamiltonian of the system can be written in the 
rotating wave approximation 6 as following: 
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Sb^bk + AB^Bk + T k ( BlK + B k b 



(2) 



Here 8 = Ejjc mi — fiw and A = Eis — huj are the detunings 
of the laser frequency from the FE and WE respectively; 
N is the number of unit cells in the organic; E is the elec- 
tric field of the incoming light; £> k and fo k are annihilation 
operators for the WE and FE respectively. 

Because the frequency of the laser pulses is much smaller 
than u, one can consider the laser field turned on adiabat- 
ically. Then, the equations of motion for any operator can 
be solved in a very straightforward manner. Apply to that 
operator the unitary transformation that diagonalizes (2) 
and then introduce an imaginary part to all detunings, pos- 
itive for creation operators and negative for annihilation 
operators. Since the Hamiltonian involves only Bose op- 
erators and is quadratic, the unitary transformations that 
diagonalizes it can be obtained in an analytic form. The 
diagonalization procedure can be written as: 
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° To assure thermalization of the hybrid excitons we propose to 
pump the sample with consequent laser pulses coming one after 
another withing the hybrid life-time. 

6 The rotating-wave approximation has been invoked by eliminating 
the anti resonant terms of the interaction with light; this approx- 
imation limits the calculation only to the resonant features of the 
response. 



where Si = (x^bl - x^b^j + (yk-B k - y^B^j eliminates 

the linear terms in + b^j . The transformation coeffi- 
cients age given by 



AS-Tl 



AS-rl 



The bilinear cross term B^bk + B^b^ can be diagonalized 
away through a rotation by (/> k in the FE-WE coordinate 
space [Mukamel (1995)] of the form: 

S 2 = 0kBk& k - faBlbk, tan (20k) = — ^ 
The overall diagonalized Hamiltonian has the form: 



4r k } BtB k + 




(3) 



The factors in front of number operators B^B^ (6 k 6 k ) 
represent the upper branch and lower branch of the hy- 
brid at exact resonance between FE, WE and the photon 
Eis = Edcm2 — The dispersion is illustrated on Fig. 2 
[Roslyak and Birman ( arXiv:cond-mat/0703268)]. 
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Fig. 2. The solid lines represent upper and lower branches of the 
quadrupolc-dipolc hybrid dispersion when the coupling is calculated 
in the parabolic approximation and the induced Stark effect is taken 
into account; the dash lines correspond to the coupling parameter for 
different approximations: (1) - infinite IQW and the induced Stark 
effect is neglected, (2) - parabolic approximation, (3) - infinite IQW 
with the induced Stark effect treated as perturbation 

With the use of this unitary transformation, the direct 
contribution of the exciton transition to the induced polar- 
ization per unit area of the interface is 

(4) 



Q 



Here the initial state \i) corresponds to the FE exciton, as 
we neglected the quadrupolc-light interaction. The hybrid 
life-time is introduced into the expression above through 
the complex part of the detunings S — > 5—ih>y, A — > A—ih-y 
as discussed above. 

Also, due to the localization of the Frenkel non-resonant 
excitation, we take the FE into account as corrections to 
the local field. Assuming the organic to be an isotropic 
medium, one can write total electric field to be 
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The second term represents the local field contribution, 
where a non - res is the background polarizability (defined 
by non-resonant excitations in the organic) , while a res is 
the direct exciton contribution to the polarizability corre- 
sponding to ( 4 ). The last equation can be rewritten as: 
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1 - (4™ res /3) (e + 2) /3 
where the Lorentz-Lorenz relationship 
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has been used. Thus, the local field corrections results in 



the exciton frequency shift 5 
susceptibility is given by: 
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Taking into account that N/ S : 
efficient is given by: 



1/ap, the absorption co- 
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Let us now make an important observation, that this ab- 
sorption coefficient is fundamentally different from that for 
the case of the DDH. In case of the QDH, one has TvyF 3> 
ft-lhyb 3> hjiS while for a dipolc-dipolc hybridization they 
are of the same order of magnitude. So, when A — > the 
absorption for the QDH is zero even for a smallest coupling. 

For the QDH using the given absorption coefficient (5), 
the transmission through the sample of an arbitrary thick- 
ness (See Fig.l) can be written as: 



Tr = 1 — exp 



(jhAY 
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I 



(6) 



Here z' can be viewed as the width of the narrow strip 
around the interface where the hybridization occurs and 
I = "^f^F" is the hybrid absorption length. 



Due to the big oscillator strength and narrow line width 
of the hybrid the absorption length I is extremely small. 
This is balanced out by the fact that the hybridization oc- 
curs only in a narrow strip around the interface between 
the organic and inorganic. The rest of the crystal is approx- 
imately transparent. So we can assume that the absorbing 



region is equal to the Bohr radius of the hybrid, i.e. z' = 
as = 7 A, therefore z'/l w 1. The value of the wave vec- 
tor k is controlled by the oblique angle of incidence for the 
probe laser. 

The "weak" interaction of the WE in the quantum well 
with the LA phonons leads to a red shift of the semicon- 
ductor energy gap and therefore to the red shift of the WE 
band gap [Passler (1998); Weiser (1992); Varshi (1967)] 
and does not affect the binding energy of the IS* exciton 
[Snoke et al. (1990)]. Hence, for small deviation of the ex- 
citon gas temperature from the helium bath temperature 
one can consider the shift to be linear [Roslyak and Birman 
(arXiv:cond-mat/0703650v2); Dasbach (2004)]: 



E w (T)=E 1S (0)-K(T-T bath ) 



(7) 



K 



where k = 0.3 

The final temperature change given by the laser heating 
and helium cooling can be written as: 

rlT 

— = HP ex [l-Tr( T ,T)}-C[T-T bath ] (8) 

ClT 

Here the phenomcnologically introduced constants H = 



5.0 x 10 6 Ks- x W- x and C = 2.8 x 10 3 



giving the heat- 



ing and cooling rate, respectively [Dasbach (2004)]. The 
transmission (6) depends on the temperature through the 
detuning A (T) according to the equation (7) . The equi- 
librium temperature is given by the stationary solutions of 
the equation (8). This will govern the temperature and ab- 
sorption bi-stability. 

3. Results and discussion 

The QDH structure offers a more complicated picture 
of the nonlinear photo-thermal effect then just a slab of 
cuprous oxide itself. As we now show it is not a mere en- 
hancement of the bi-stability effect but the appearance of 
two or more well pronounced bi-stability regions in tem- 
perature and transmission (absorption) spectra. In case of 
the QDH the photo-thermal rise of the temperature leads 
not only to the red shift of the IS exciton energy band but 
to an effective breaking of the resonance with the FE and 
reduction of the hybridization. 

Many properties of the temperature, absorption and ef- 
fective dispersion of the hybrid system can be obtained 
even before a numerical solution for the system of equa- 
tions (7,8). The right hand side of the equation (8) can be 
visualized as an implicit function of the temperature and 
detuning (See Fig. 3). The actual temperature of the system 
depends on the way one changes the frequency of the probe 
laser. If one lowers the frequency starting from the detun- 
ing 2 meV (point L on Fig. 3) off resonance, then the sys- 
tem follows the path i?^1^3^4^5^0^7^L 
to assure minimum temperature change with the changing 
detuning. If one increases the frequency starting from the 
detuning —2 meV (point R on the Fig. 3), then the system 
follows the path i^7^6^5^4^2^1^i?. 
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Fig. 3. The temperature of the sample in K, versus the detuning in 
units of the hybrid exciton line-width *yh. The coupling parameter 
is taken to be equal to the line width of the hybrid and the subplot 
(b) corresponds to the Stark reduced coupling. The thick solid curve 
corresponds to the laser intensity P ex = 3 mWt and maximum 
temperature 5 K, and another thin curve stands for the intensity 
P ex = 1 mWt and maximum temperature 2.5 K 

On the other hand one can consider bi-stability as a pos- 
sibility for the the system (7,8) to have more then one sta- 
tionary solution. This non-linearity grows with P ex . To il- 
lustrate the concept it is convenient to plot the heating and 
cooling rate as a function of the temperature. Then the sta- 
tionary temperature is given as the intercept points ( See 
Fig.4). 
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Fig. 4. The numerical solution of the nonlinear equation (8) is given 
crossing of the cooling (straight line starting from the origin) and 
the detuning dependent photo-heating. The figure a)-k) correspond 
to the different laser detuning from the low temperature limit of the 
15 quadrupole exciton as ^- = 5; 3.7; 3; 2.5; 1.8; 0.5; 0.4; 0; -0.8 
in units of the damping parameter. 



Now let us focus on the different branches of the hybrid. 
It is clear from the Fig. 3 and Fig. 5 that the response from 
two branches is asymmetrical. Further on we are going to 
refer to the lower branch as Wannier-like and to another 
as Frcnkcl-like. Both of them are defined by the interplay 



between the thermo-induced red shift of the WE exciton 
k (T — Tb a th) and the excitation detuning A. 

Indeed, the maximum of the absorption given by the 
poles of the expression (5): 
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WE and FE is dominant: 



k (T — Tbath) <C 4r^, then both branches are symmet- 



rical and the maximum of absorption corresponds to the 
usual hybrid dispersion S± = ±IV Otherwise, strong 
asymmetry reveals itself. The case <5+ — > k (T — Tbath) de- 
fines the maximum of Wannier-like bi-stability, dominated 
by the temperature induced red shift. The Frcnkcl-like bi- 
stability corresponds to 5- — > and is dominated by the 
hybridization effect. 

For the up-scanned Wannier-like branch when the exci- 
tation energy 6 catches up with the escaping exciton res- 
onance S + the actual temperature cut-off (Path 2 — ► 1 
on the Fig. 3 and Fig.4, (1)) occurs when the laser heating 
reaches its maximum <5 + = A. When the excitation energy 
increases further the absorption reduces and the sample 
cools rapidly. On the down-scan the sample is cold even 
beyond the cut-off energy (the temperature is close to the 
lower cross-point on Fig.4, (1)). 

In the same fashion the temperature rises rapidly on 
the Frenkcl-like branch (Path — > 7 on the Fig. 3 and 
Fig.4, (b)) with a down-scan. The temperature goes down 
and the detuning between FE and WE reaches its mini- 
mum and an abrupt heating occurs. For the up-scan the 
sample is heated enough (close to the upper cross-point on 
Fig.4, (b)) to maintain large detuning between the FE and 
WE but drops down rapidly when the FE and WE detun- 
ing reaches its maximum <5_ — A (Sec Fig.4,(d)). 

The bi-stability effect and the temperature rise of the hy- 
brid is greatly enhanced (meV) compared to bulk cuprous 
oxide itself (peV) mainly due to the fact of gaining oscilla- 
tor strength from the organic part of the hybrid. 

Numerical search for a stationary solution of equation 
(8) reveals a fine structure of the multiple bi-stability (See 
Fig. 5, Fig. 6). For the up-scan regime the distance between 
the solution branches 4^2 and 3 — > 2 on the Fig. 3 
can become smaller then n (T — Tbath) and the temperature 
drops down on the branch 3^2. For the same reason when 
the intensity is big enough and the coupling between the 
FE and WE is reduced by the Stark effect fine structure 
can appear on the branches 5 — > and 6 — > on the down- 
scan. 

The model described above is valid for small tempera- 
ture raise from the Tbath only. But it can be generalized for 
bigger temperature change by using a more elaborated non- 
linear temperature dependence of the WE exciton [Varshi 
(1967)] and possible effect of the optical phonons [Passler 
(1998)]. Although, when one increases P ex to the point 
when the density of the hybrid is close to the saturation 
density then other nonlinear effects due to exciton-exciton 
scattering and possible condensation start to play a con- 
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Fig. 6. The temperature of the interface between organic DCM2:CA:PS and inorganic cuprous oxide QW for different detuning of the laser 
field from the energy of the 15 quadrupole exciton take at T = Tbath and different wave vectors k of the hybrid. The blue color corresponds 
to the temperature of the surrounding helium 1.7 K. The graphs represents the up (A) and down (B) laser scanning. The stark effect is taken 
into account on the up-scan (C) and down-scan (D) graphs. 
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Fig. 5. The numerical results for the temperature of the sample in K, 
versus the detuning in meV, P ex = 3 mWt. The coupling parameter 
is taken to be equal to the line width of the hybrid. The bi-stability 
is a result of multiple solutions and reveal itself as a sudden drop or 
raise in the temperature depending on the direction of the scan. 

siderable role. Such non-linear effects are presently under 
investigation. 
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